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ABSTRACT 

This  thesis  attempts  to  develop  a  program  for  the  Whirlwind 
computer  operating  in  the  interpretive  mode  for  the  flux  distribution 
in  a  finite  cylindrical  reactor© 

The  reactor  chosen  for  investigation  was  a  hypothetical  organic 
moderated^,  highly  enriched^  cylindrical  reactor  having  a  power  output 
eapable  of  fulfilling  one  half  the  propulsive  requirements  of  a 
nuelear  powered  naval  cruiser o 

H«at  transfer  and  basic  physics  calculations  were  accomplished 
to  determine  preliminary  size  and  material  requirement  of  such  a 
r«aetoi>o 

The  flux  equations  were  derived  using  numerical  integration 
methods  suitable  for  computer  programming* 
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Conclusion 

1)  The  code  as  presently  programmed  is  incorrect* 

2)  The  grid  spacing  as  used  is  much  too  coarse.  This 
difficulty  is  not  inherent  in  the  coded  program* 

3)  The  mode  of  computer  operation  is  unsuitable  for  an 
iterative  process  such  as  this  because  it  results  in 

■  a  code  that  progresses  too  slowly© 

Recommendations 

1)  Prior  to  the  use  of  any  computer  a  formal  or  semi-formal 
©ourse  in  programming  should  be  completed© 

2)  A  grid  spacing  comparable  to  the  mean  free  path  of  the 
neutron  should  be  selected© 

3)  The  fastest  possible  mode  of  calculation  on  the  computer 
assailable  should  be  usedo 

k)  For  a  eoded  program  of  this  difficulty  a  faster  machine 
with  a  much  larger  fast  memory  is  very  desirable© 


Thesis  Supervises1  Melville  Clark.  Jr. 

Title  Assistant  Professor  of  Nuclear  Engineering 
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INTRODUCTION 

The  desirability  of  a  uniform  power  distribution  (P   /P  )  ■  1 

max  av 

has  long  been  realized  since  it  will  not  only  decrease  the  mass  of 
the  fuel  required  for  the  attainment  of  criticality  but  also  will 
increase  the  total  power  output  by  allowing  the  whole  core  instead 
of  only  a  small  fraction  of  the  reactor  to  operate  at  the  highest  allow- 
able temperature*  Operation  at  the  highest  possible  temperature  allowed 
by  metallurgical  considerations  over  the  entire  reactor  will  increase 
the  average  core  life  since  more  uniform  burn  up  is  attained* 

It  is  possible  to  achieve  this  uniform  power  distribution  by 
various  means •  One  of  the  first  methods  to  be  used  was  that  of 
poisons*  In  this  method  long  lasting  poisons  are  placed  in  the 
reactor  structure  in  the  anticipated  locations  of  high  flux*  By 
placing  poison  in  these  regions  the  peak  flux  is  reduced  to  such  a 
level  that  the  power  density  is  approximately  uniform*  This  pro- 
ca&Ha?e=is  a  very  wasteful  one  from  the  point  of  view  of  neutron  economy. 
The  regions  of  the  reactor  that  are  contributing  greatly  to  the  overall 
neutron  efficiency  of  the  reactor  are  penalized  to  the  extent  that 
they  can  be  no  more  efficient  than  the  most  inefficient  region*  For 
power  producing  reactors,  particularly  mobile  ones,  this  waste  of 
neutrons  is  very  undesirable* 

The  use  of  poison  is  the  first  of  two  methods  that  exploit  a 
nonuniform  macroscopic  cross  section  to  achieve  the  uniform 
distribution  of  power*  The  second  consists  in  the  use  of  a  nonuniform 
macroscopic  fissiqn  cross  section*  This  method  allows  each  region  of 


the  reactor  to  operate  most  efficiently.  The  efficiency  is  determined 
by  the  position  of  the  region  within  the  reactor©  The  leakage  from 
a  reactor  causes  a  reduction  in  the  value  of  the  flux  at  large 
distances  from  the  center  of  the  core*  To  compensate  for  this  reduc- 
tion in  flux  the  macroscopic  fission  cross  section  is  increased  in 
the  Regions  far  from  the  center  and  decreased  in  the  central  portions. 
To  this  endc,  the  fuel  concentration  is  increased  near  the  edges* 

Neutron  economy  can  also  be  improved  by  surrounding  the  core  with 
neutron  reflector  material.  A  light  substance  having  a  relatively 
low  absorption  cross  section  performs  this  function  very  well5  i.e., 
heavy  water  or  beryllium. 

Since  the  reactor  to  be  investigated  is  intended  for  a  hypothetical 
naval  cruiser 5  the  use  of  poisons  to  achieve  a  uniform  power  distribu- 
tion is  not  desirable,,  for  a  compact  reae  tor  is  required.  Though  the 
reactor  must  be  compact,,  it  is  desirable  to  install  a  reasonably  thick 
reflector  for  the  purpose  of  improving  the  neutron  economy  s  attenuating 
^©^eutrons  leaking  out5  and  cooling  of  the  reactor  shell  and 
pressure  vessels 

It  would  be  possible  to  have  a  nonuniform  macroscopic  cross 
section  of  the  reactor  in  each  of  its  directions*  For  the  spherical 
reactor  the  variation  would  be  radial  only5  but  the  application  of  a 
spherical  reactor  to  shipboard  use  is  impractical.  The  norfyal  con- 
figuration of  spaces  aboard  a  ship  is  rectangular  and  a  spherical 
reactor  does  not  adapt  very  well  to  this  arrangement.  However  <,  a 
rectangular  reactor  allows  efficient  utilization  of  space  aboard  a 
shipj,  but  from  the  nuclear  point  of  view  this  shape  is  not  too  desir- 


able  because  leakage  depends  on  the  relationship  of  surface  area  to 
volume.  It  can  be  seen  that  the  neutron  economy  decreases  in 
progression  from  the  sphere  to  the  cylinder  to  the  rectangular 
parallepiped© 

With  these  considerations  in  mind  a  right  cylindrical  reactor 
shape  will  be  investigated*  With  this  shape  it  would  be  possible 
to  vary  the  properties  in  two  directions ,  radially  and  axially© 
While  this  goal  is  theoretically  interesting,  it  is  prohibited 
when  cost  is  considered*  The  reactor  that  is  to"be  investigated  will  have  a 
uniform  composition  in  the  axial  direction  with  finite  radial  segments 
differing  in  properties© 

For  the  shipboard  use  of  a  nuclear  reactor  power  plant  it  is 
desirable  to  have  a  plant  as  light  and  safe  as  can  possibly  be 
designed  within  the  space  allowed  but  still  without  incorporating 
undue  quantities  of  the  "exotic"  and  expensive  materials  or  coolants. 
One  possible  solution  is  the  use  of  a  highly  enriched  reactor  using 
a  coolant  that  permits  high  operating  temperatures  at  low  pressures 
and  does  not  become  dangerously  radioactive  while  passing  thru 
the  reactor  core© 

A  coolant  that  has  a  low  vapor  pressure  at  high  operating 
temperature  will  allow  a  system  pressure  which  is  much  lower  thah  that 
employed  in  the  present  day  water  cooled  and  moderated  shipboard 
reactor  installation*  The  use  of  such  a  coolant  will  not  only  result 
in  much  better  steam  conditions,  but  it  will  also  make  possible  an 
overall  weight  reduction,  since  a  low  pressure  system  can  be  used© 
The  saving  in  weight  is  particularly  significant  in  the  pressure 
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■vessel© 

If  the  coolant  does  not  become  radioactive  in\',its  passage  thru 
the  reactor  corej,  the  shielding  can  be  concentrated  primarily  around 
the  reactor  proper*  provided  that  impurities  in  the  coolant  are 
kept  at  a  very  low  concentration,,  Even  with  a  small  amount  of 
impurities,,  the  radiation  emitted  by  the  coolant  would  be  low 
enough  to  permit  reduced  thicknesses  of  shielding  around  the  primary 
loop* 

A  liquid  organic  hydrocarbon  will  fulfill  the  requirements  listed 
above s  but  also  will  present  some  unsolved  problems.  It  will  (l)  tend 
to  increase  the  size  of  the  reactor  by  decreasing  the  hydrogen  modera^ 
ti on  relative  to  a  water  cooled  and  moderated  reactor.,  (2)  tend  to 
increase  the  rea@tor  size  due  to  its  less  efficient  heat  transfer 
capabilities  compared  with  water  and  (3)  exhibit  both  thermal  and 
radiolyti®  breakdown  in  which  the  original  hydrocarbon  will,,  primarily 
by  polymerization*)  change  into  longer  chain  hydrocarbons  with  the 
evolution  of  relatively  large  amounts  of  gas  and  a  great  increase  in 
viscosity  of  the  coolant© 

The  latter  problems  will  be  investigated  at  the  National  Reactor 
Testing  Station  ©n  the  Qrganic-Moderated  Reactor  Experiment  (ORME)© 
This  experiment  is  not  a  prototype  for  an  organic  power  reactor  but 
is  being  built  to  test  the  feasibility  of  its  use  as  a  reactor  coolanto 
At  the  West  Milton  Annex  of  the  Knolls  Atomic  Power  Laboratory  the 
Atomic  Energy  Commission  plans  to  investigate  the  feasibility  of  an 
organic  moderated  and  cooled  reactor  for  naval  propulsion.  This  will 
be  conducted  on  the  Naval  Organic  Reactor  Experiment s   (NORE)© 
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The  heat  transfer  properties  of  the  various  liquid  organic 
hydrocarbons  that  are  presently  under  consideration  are  not  generally 
available*  Some  presently  under  investigation  are  diphenyl, 
0-terphenyls  M-terphenyl  and  P-terphenyl.  A  proprietary  heat  transfer 
medium  "Dowtherm"  is  available  commercially  and  is  composed  of  73«£# 
diphenyl  oxide  and  26»5%  diphenyl.  It  has  a  working  range  of  2i£0- 
75>0 °F«  at  pressures  less  than  lU5  psig.  Since  "Dowtherm"  closely 

resembles  the  hydrocarbons  that  will  be  investigated  on  the  QRME  and 

[2] 

its  heat  transfer  properties  are  readily  available   it  will  be  used 

in  this  investigation. 

The  fuel  element  chosen  for  this  investigation  is  of  the  MTR 

[3l 

(Materials  Testing  Reactor)  type.  "   An  effective  length  of  element 

was  set  at  li8",   his  length  was  chosen  so  that  the  overall  height 
of  the  reactor.,  control  devices  and  shield  would  fit  into  the  present 
day  engineering  spaces  of  a  cruiser.  Any  additional  height  would  have 
an  adverse  effect  on  stability,  since  it  would  result  in  a  higher 
shield  deck  and  thus  an  increase  in  the  height  of  the  center  of 
gravity  of  the  ship. 


Fuel  Element  Dimensions 
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PROCEDURE 
Heat  Transfer" 

The  controlling  variables  in  this  study  were  (l)  the  maximum 
surface  temperature  of  the  fuel  element  or  the  maximum  temperature 
at  the  center line  of  the  fuel  element  and  (2)  the  maximum  coolant 
speed.  The  maximum  surface  temperature  of  the  fuel  element  may  be 

controlling,,  since  the  preliminary  investigations  on  the  ORME  have  .   . 

[h] 

shown  that  above  800°  F  the  thermal  breakdown  is  increased  greatly. 

However  9  liquid  loading  requirements  of  a  naval  vessel  make  storage 
of  used  coolant  a  simple  matter  and  in  the  interest  of  higher 
efficiency  some  breakdown  will  be  accepted.  Therefore,  a  fuel  element 
surface  temperature  less  than  8^0°  F  will  be  accepted.  If  the  maximum 
fuel  temperature  at  the  centerline  is  allowed  to  exceed  a  value  of  about 
1130°  F5  a  phase  transformation  of  uranium  is  encountered  in  which  an 
anisotropic  expansion  and  possible  rupture  of  the  fuel  element  may 
result©  Therefore «,  a  second  design  criterion  will  be  required  that 
the  fuel  element  centerline  temperature  be  less  than  1130°  F. 

Preliminary  studies  on  the  ORME  have  also  shown  that  coolant 
speeds  greater  than  l£  feet  per  second  result  in  excessive  erosion. 
As  a  third  design  criterion^  it  will  be  required  that  coolant  speed 
be  less  than  l£  feet  per  second. 

The  first  investigation  was  with  a  reactor  coolant  outlet  temper- 
ature of  75>0°  F  with  various  temperature  rises  across  the  reactor. 
This  investigation  was  amplified  by  varying  coolant  speed  while  holding 
the  temperature  rise  across  the  reactor  constant.  Each  of  these 
situations  were  further  studied  by  considering  the  MTR  fuel  element 
modified  to  the  extent  that  the  plate  spacing  was  narrowed  to  80$ 
and  60$. 

8  Nomenclature  for  Heat  Transfer  on  pages  22  and  23. 
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Though  it  was  known  in  advance  that  the  speed  would  be  kept  at  or 
less  than  i£  feet  per  second  calculations  were  pursued  with  higher 
values 0  This  was  done  to  show  that  a  reduction  in  size  can  be  accom- 
plished in  the  future^  if  the  erosion  problem  at  high  coolant  speeds 
can  be  overcome 0 

A  heat  rate  requirement  of  6«06  x  10  Btu/hr  was  chosen  based  on 
the  full  power  steaming  rate  of  a  typical  cruiser* 

With  this  value  of  Q  (6«06  x  10  )  a  coolant  outlet  temperature 9 
coolant  temperature  rise5  coolant  speed  and  plate  spacing  are  chosen* 
Since  the  inlet  and  outlet  temperatures  are  known  the  average  values 
of  density  and  specific  heat  may  be  calculated*  Formulae  19  2S   3$  an^ 
k  are  used  respectively  to  calculate  w5  A  9   G  and  n. 


w  *  -s-2—  .  Q  A  (1) 

C  AT      f 
P 


w 


f   3600fv 


(2) 


G-^-5cooT  (3) 

f 
n  «  —  (k) 

Since  the  dimensions  of  the  '^uel  element  for  each  case  are  known, 
both  the  hydraulic  diameter  and  the  total  flow  area  can  b)f  calculated© 
With  the  product  of  the  mass  flow  rate  and  the  hydraulic  diameter  th© 
produce  of  the  heat  transfer  coefficient  and  the  hydraulic  diameter  can 
be  read  from  the  graph  in  the  "Dowtherm"  booklet.  Formulae  55  6,,  7$ 
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Datum  from  the  calculation  of  a  coolant  outlet  temperature  of  750°  F 
for  various  values  of  temperature  rise  and  speed  for  cases  1^  2  and  3 
are  in  Tables  I5  11^  and  HI  and  Figures  19   II,  and  HI  respectively* 
The  cladding  temperature  existing  for  each  of  these  conditions  is 
plotted  on  Figure  IV, 

The  narrowest  plate  spacing^,  case  3$   gives  a  significantly  smaller 
reactoro  Since  the  pumping  power  of  any  of  these  cases  proved  to  be 
very  small  compared  with  the  total  power  output^  it  was  decided  to 
investigate  the  result  of  a  decrease  of  coolant  outlet  temperature  fo? 
this  case  only*  Coolant  outlet  temperature  was  reduced  to  700°  V  and  the 
speed  was  held  at  1$   feet  per  second©  Datum  from  this  calculation  are 
in  Table  IV©  These  data  were  plotted  as  a  dotted  line  on  Figure  IV  • 

For  these  twelve  conditions  the  fraction  of  uranium  in  the  uranium^ 
zirconium  dispersed  fuel  plate  was  calculated  and  plotted  on  Figure  V« 
It  is  desirable  that  the  fuel  be  of  such  a  composition  that  it  contributes 
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TABLE  I 

Datum  for  coolant  outlet  of  750°  F 
and  wide  plaifce  spacing 
(case  1) 

A.  REACTOR  AT  -  200  °F 

v         G   2         h  2      Af         n 
(ft/sec)  (lb/see-ft  )  (Btu/hr-ft  -°F)  (ft2)  (no.  of  plates) 

5  236,1  636.  5.36  29685 

15  709.0  1611.  1.79  896 

25  1180.0  231*0.  1.07  538 

50  2360.0  1*090.  0.51*  269 

v  0^  x  108         T.-T-,^     TrT2     T2-Tb       T2  fu 

(ft/sec)     (Btu/nr-ft  )     (eF)       (°F)       (°F)       (°F)  % 

15  1**86  20.7      62.2      215.8    1015.8*      25.1* 

B.  REACTOR  AT  =  150°F 


(ft/sec) 

G 

2 
(lb/see-ft  ) 

h 
2 

(BtuAr»ft  - 

Af 
•F)  (ft2)  (no. 

n 
of  plates) 

5 

231.9 

682. 

7-26 

35639 

15 

705.0 

1591. 

2.1*2 

19211 

25 

1160.0 

2320. 

1.1*5 

727 

50 

2318.0 

1*01*0. 

0.73 

361* 

(ft/sec) 

Qv  *  108 
(BWhr-ft  ) 

T°"T1  T1~T2 
(°F)   (*F) 

VTb   T2 
(°F)   (*F) 

% 

15 

3.60 

15.1*   1*6.0 

197.3  91*7.3* 

19.5 

excessive  temperature 

Y  1 

'uranium  content  too  low 


Table  I  (Continued) 


C.  REACTOR  AT  -  100  °F 

v  G                             h  Af  n 

(ft/sec)  (ib/see-ft  )     (Btu/hr-ft  -°F)     (ft2)     (no.  of  plates) 

5  227o                         669.  11.01  55520 

15  681*                       1569.  3.67  1,839 

25  1135.                         2295.  2.02  1,103 

50  2270                           1*000.  1.10  552 

▼  0^  x  vfi         To-T-t     TrT2  T2«Tb       T?  fu 

(ft/sec)  (Btu/hr-ft  )     (°F)       (eF)  (°F)       (°F)  % 

15  2.37            10.1      30.3  133.8    883.8*  13.5 

D.  REACTOR  AT  -  50  eF 

v  G                             h  Af  n 

(ft/see)  (Ib/see-ft  )     (Btu/hr-ft  -*F)     (ft2)     (no.  of  plates) 

5  221.8                           6#.  22.1*0             11^230 

15  667.0                     152*5.  T.\6  3,730 

25  1110.0                        2275.  2*.U7  2,21*0 

50  2218.0                       3910.  2.21*  1,123 

v  Ox  10           T  =T       T  -T  T  «=T         T  f 

(ft/sec)  (Btu/hr~ft  )     (°F)       (®F)  (°F)       (°F)  Jg 

15  1.17           li.97      15.0  66.2      816.2  6.99 


excessive  temperature 

Y 

uranium  content  too  low 
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TABLE  H 

Datum  for  coolant  outlet  of  7£0°F, 
coolant  speed  of  15>  ft/sec,  and  medium  plate  spacing 

(case  2) 

AT  G  h  n  Qv  x  108 

(*F)     (lb/sec-ft  )     (Btu/hr-ft  -°F)     (no.  of  plates)     (Btu/hr-ft3) 


200. 

709. 

1671. 

1121. 

3.89 

150. 

705. 

1660. 

1515. 

2.88 

100„ 

681. 

1639. 

2300* 

1.89 

50. 

667. 

1631. 

1*660. 

0.91* 

AT 

T,-^ 

T  <=T 
1  2 

T  -T 
2  b 

T2 

*u 

(eF) 

(°P) 

(•F) 

(°F) 

(*F) 

% 

200 

16«6 

ltf.7 

193.9 

9l43-9* 

19.6 

150 

12.3 

36.9 

lhk*2 

891*.  0 

16.1 

100 

8.1 

2l*.3 

96.3 

81*6.3 

11.0 

50 

1*.0 

12.0 

1*7.6 

797.6 

5.6Y 

* 

temperature  excessive 

Y 
uranium 

content 

too 

low 
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TABLE  ni 

Datum  for  coolant  outlet  of  750°F, 
coolant  speed  of  15>  ft/sec^,  and  narrow  plate  spacing 

(case  3) 

AT        G  h  n  n         Q  x  108 

2  2  ^3 

(*F)  (lb/sec-ft  )  (Btu/hr-ft  -°F)  (no*  of  plates)  (Btu/hr-ft  ) 


200 

709 

1796 

lii96 

2.9 

150 

705 

1775 

2020 

2.2 

100 

681 

1731 

3061 

Uk 

50 

667 

1718 

6230 

0.70 

AT 
(°F) 

(«F) 

T  ~T 
(•F) 

T2 

(ep) 

f 

200 

12.1|. 

37.3 

135.2 

885.2* 

16.3 

150 

9.2 

27.6 

I01e3 

851.3 

120ii 

100 

6el 

18.2 

68.5 

818.5 

QokX 

50 

3o0 

9o0 

33.9 

783.9 

1*.3Y 

* 
temperature  excessive 

Y 
uranium 

content 

too 

low 
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Volant  sp*«£  is,  ft/soe 
13 


nu^Obtc  of  plat-as 


Coolant  speed  in 
Hi 


of  plates 


number  oi  pl»t*g  -J 


.id. 


TABLE  IV.  . 

Datum  for  coolant  outlet  of  700 *F, 
coolant  speed  of  15  ft/sec*  and  narrow  plate  spacing 

(case  3) 

AT        G  h  n         Qy  x  108 

(«F)  (lb/sec-ft  )  (Btu/hr-ft  -°F)  (no.  of  plates)  (Btu/hr-ft3) 


200 

739.5 

1826 

871 

5.0 

150 

726„8 

1812 

1161* 

3.7 

100 

712.5 

1775 

1767 

2.5 

AT 
(eF) 

VT2 

T  -T 
1  2 

(•F) 

T  -T 
(*F) 

(*F) 

fu 
* 

200 

21.3 

61^0 

228.1^ 

92U.1; 

26.0 

150 

15.9 

1*7.9 

172.2 

872.2 

20.3 

100 

10.5 

31.5 

115.8 

815.8 

11.0 

temperature  excessive 

Y 
uranium 

content  too  low 
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W^ 


to  the  structural  strength  of  the  reactor.  A  uranium  content  of  10J6  by 
weight  or  greater  has  desirable  strength  characteristics. 

Of  the  six  conditions  that  have  a  fuel  element  surface  temperature 
of  less  than  8£0°  F  only  two  have  a  uranium  fraction  of  above  10#. 
These  are  shown  in  Table  V. 


t 

CO 

T 

case 

n 

T2 

750 

100 

3 

3061 

818 

700 

100 

3 

1767 

81^ 

Table  V  Particularly  Interesting  Cases 

Figure  "VI  shows  the  variation  in  number  of  plates  for  case  3 
holding  the  coolant  speed  at  15  feet  per  second  and  varying  the 
temperature  rise  across  the  reactor. 

The  size  resulting  from  the  condition  with  a  coolant  outlet 
temperature  of  700°  F  is  better  for  the  following  reasonst 

1)  it  is  by  far  the  smaller, 

2)  the  increase  of  uranium  fraction  by  2  1/2^  gives  a 
significant  increase  in  the  strength  of  the  fuel  region. 

With  this  condition  the  fuel  element  surface  temperature  will  be 
less  than  the  prescribed  value  of  850°  F.  The  amount  of  damage  that  will 
result  to  the  coolant  from  this  temperature  and  the  accompanying 
radiation  is  uncertain.  However,  for  shipboard  applications  it  does 
not  appear  to  be  a  serious  penalty. 

In  the  present  day  cruiser  design  for  torpedo  protection  there  is 
available  a  large  amount  of  bunker  fuel  oil  to  partially  fill  the 
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layers  of  the  torpedo  protective  system.  Certain  layers  of  the  torpedo 
protective  system  are  required  to  be  liquid  loaded.  This  system  could 
be  used  for  the  storage  of  the  organic  moderator  coolant.  The  unradiated 
coolant  could  be  stored  in  the  inner  layers  of  the  system  and  that 
which  had  been  irradiated  would  be  stored  in  the  outer  layers  to  protect 
the  crew  from  any  induced  radioactivity.  The  long  chain  hydrocarbons 
would  tend  to  settle  to  the  lower  parts  of  the  tanks  where  they  could 
be  removed  during  periods  of  overhaul.  With  the  installation  of  a 
purification  system  this  process  of  separation  might  be  handled  con- 
tinuously aboard  ship  in  the  same  manner  that  lubrication  oil  is 
purified  and  recycled. 
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Nomenclature 


2 
A        flow  total  area,  (ft  ) 


f 

A  flow  area  per  plate  case  n,  (ft/plate) 

A  heat  transfer  area,  (ft/plate) 
h 

c  average  specific  heat.  (Btu/lb  °F) 
P 

D1  hydraulic  diameter,  (inches) 

f  fraction  of  uranium  in  fuel  region 
u 

G  mass  flux,  (lb/sec-ft  ) 

h  film  coefficient,  (Btu/hr-ft  -°F) 

k  thermal  conductivity  of  cladding.  ( Btu/hr-ft- °F) 
c 

k  thermal  conductivity  of  cladding,  (Btu/hr-ft-*F) 
f 

L  fuel  length,  (inches) 

n  number  of  fuel  plates 

3 

Q  total  heat  generation  rate,  (Btu/ft  ) 

3 

Q  heat  generation  rate  per  unit  volume  of  fuel,   (Btu/hr-ft  ) 

t  fuel  cladding  thickness,  (inches) 
c 

t  coolant  outlet  temperature,  (*F) 

CO 


t        fuel  thickness,  (inches) 


m 


t        spacing  between  plates,  (inches) 


s 


T        bulk  temperature  of  coolant,  (*F) 
b 

22 


Nomenclature  (Continued) 

To  fuel  center line  temperature,  (®F) 

T  fuel-clad  interface  temperature,  (°F) 

T  fuel  element  surface  temperature,  (*F) 

U  overall  heat  transfer  coefficient,  (Btu/hr-ft  -T) 

v  coolant  speed,  (ft/sec) 

w  mass  rate  of  flow,  (lb/hr) 

W  fuel  width,  (inches) 

AT       overall  temperature  difference  across  reactor,  (*F) 

3 
p       average  density,  (lb/ft  ) 

©       film  temperature  drop,  (*F) 
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Basic  Reactor  Physics 

Criticality  calculations  were  made  using  a  on*  group  reflected 
homogenous  model  with  case  3>  narrow  fuel  elements.  Based  on  heat 
transfer  considerations,  inlet  and  outlet  coolant  temperature  were 
set  at  600°  and  700°  F,  respectively,  with  a  coolant  speed  of  15>  ft/sec. 
The  assumptions  in  this  calculation  were 

1*  The  structural  members  inside  the  reactor  were  of 
type  3li7  stainless  steel  and  composed  ten  percent  of 
the  total  reactor  weight, 

2.  A  fast  fission  factor  of  unity, 

^    23£ 

3»  The  uranium  of  fuel  is  enriched  to  90%  in  U   , 

h*   A  right  circular  cylindrical  flux  distribution  of  a  cosine 

axially  and  a  Bessel  function  radially, 
5.  A  beryllium  reflector  of  236  cm  placed  around  the  reactor* 
The  reactor  was  treated  as  a  homogeneous  medium.  By  a  one  group 
calculation  the  volumes  of  each  constituent  of  the  reactor  were  cal- 
culated from  the  dimensions  given  previously  in  the  heat  transfer  sec- 
tion. These  volumes  were  then  converted  to  weights  and  the  volume 
fraction  of  the  total  reactor  as  shown  in  Table  I.  From  the  cross 
sections  given  in  BN1>32£>    effective  cross  sections  were  calculated 
for  the  operating  temperature  of  the  material.  With  these  cross 
sections  the  reactor  parameters  as  shown  in  Table  II  were  calculated. 
These  give  an  effective  multiplication  factor  of  0.89U.  This  value  is 
sufficiently  close  to  unity  that  final  adjustments  are  made  by  the 
computing  machine. 

The  primary  purpose  of  the  calculation  is  to  obtain  an  approximate 

2k 


value  of  the  critical  radius  which  is  needed  as  an  input  to  the 
computer  solution* 


TatO< 

i  I 

Reactor  Materials 

Volume 
(cm  ) 

Mass 
(Kgm) 

Volume  fraction 

Fuel 

6.95  x  1CT 

U9$ 

U 

6%$ 

0.0068 

2r 

m~.$ 

1 

y            0.1*08 

Clad  (  Zr) 

L5.01  x  1CT 

966 

Stainless  Steel 

1.99  x  1(T 

160 

Cr  GS^by  weight) 

28.Ii 

0.007 

Ni  (12%  by  weight) 

17.35 

o.ooU 

Fe  (71#by  weight) 
Coolant 

26.3  x  10 

112. 
219.5 

0.027 

[73.5  Diphenyl  oxide  by  weight  (C  H  OC.H  .)  and 

6  5    6  5 

26.5*  Diphenyl  by  weight  (M  )  ] 

C  (06.91  by  weight)  O.ltfOltf 

H  (06.03  by  weight)  0.03292 

0  (86.59  by  weight)  0.3771* 
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Tabic  II 
Reactor  Parameters 

6   «  Beryllium  shield  thickness  •     23.6     cm. 

D       -  3  £—  -  1.1761&  cm 

**  .'    i 

d      -  0.71  3  x        J?  -  33.7D  -  39.599       cm 

Atr      1-B 

2   -  2  2   V  -       0.1208  cm"1 

a    i   a.  f 

2    ©  2 

L   *'--  •      9.7387  cm 

a 

p(E)  -  «xp  %^9)lZsXw2k\  ^x  10-21*  -  0.726 

28   / 


v        ai 


a( total) 

=  /jfp  6 
R  ♦  d  /    I  H+26z 


5   a  albedo  of  Beryllium 

-B2V 
b  ka,e 

Leff  "  1+B2L2 


1.9692 

1.1*296 

0.17598 

31+3 

2 

cm 

1.332  X 

10-3 

2 

cm 

0.89I* 
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Application  of  the  Flux  Equations 

The  S  approximation  is  a  method  of  differencing  iii  which  the 
n 

independent  variable  is  quantized  into  n  intervals,  and  intermediate 
values  of  a  function  of  the  variable  are  approximated  by  a  straight 

line  between  points  at  the  boundaries  of  the  internal*  This  method 

[6] 
of  numerical  integration  was  developed  by  Bengt  G.  Carlson   • 

The  vector  flux  at  a  point  within  a  reactor  can  be  obtained  by 
the  use  of  the  S  method*  The  equations  necessary  are  derived  in 
Appendix  A,  for  two  cases  of  neutron  velocity? 

Case  A  -  The  component  of  the  vector  flux  in  the  axial  direction 
is  greater  than  its  component  in  the  radial  direction* 

Case  B  -  The  component  of  the  vector  flux  in  the  radial  direction 
is  greater  than  that  in  the  axial  direction* 

These  two  cases  are  illustrated  in  Figure  I  using  a  grid  which  is 
used  to  difference  the  transport  equation  in  the  r  and  z  coordinates, 
r  being  the  radial  coordinate  and  z  the  axial  coordinate* 

Figure  I 


The  two  cases  of  flux  at  a 


Direction  of 
neutron  motion 
Cass  B 

Direction  of 
neutron  motion 
Case  A 
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The  general  flux  equation  for  case  A  ist 

,f   '  il \  j^n ^n ^2 £ 

il      2  Vsf  h   r  / 


}t 


(1) 


where  h  is  a  value  which  is  one  greater  or  less  than  j,  depending 
upon  the  direction  of  the  radial  component  of  vector  flux,  and  n  is 
a  value  which  is  one  greater  or  one  less  than  ra,  depending  upon  the 
direction  of  the  axial  component  of  vector  flux.  Definitions  of  the 
other  quantities  can  be  found  in  the  nomenclature  part  of  this 
section. 

For  plase  B  the  flux  at  a  point  is  defined  by  the  following 
equation© 


0*m 


(2) 


where  A  for  Equation  (l)  and  (2)  is  defined  byt 


HKI-^Ww 


h  /J 


^iili^lA-^w^ 


2  \   f  h   r  J  j,ra 
h 


il      2     \ff    h      rhJ]  "h*  2  L  J,»  Mj 

^    2        _J  lj>m  &»n  h*ra  h,n__ 


(3) 
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It  should  be  pointed  out  that  the  total  macroscopic  cross 

section,  £,  is  evaluated  using  the  material  properties  at  the 

center  of  the  position  grid  of  Figure  I,  but  are  designated  by  the 

grid  values  at  the  corner  diagonally  opposite  the  point  j,», 

i.e.  j  the  point  h,n*  Further,  r  is  actually  the  mean  radius  and  is 

h 
r  *  r 
equal  to  3    h   • 

2 
The  source  at  a  point  is  defined  by  the  following  equations 


J,»   ¥    P  k*l   f  Fj,«   k-1   *  %* 


where  the  scalar  flux  0    is  equal  to 

£  H  1  P  ^/,f 

f-0  f  i-o  i    j,m 

The  fission  cross  section  and  transfer  probabilities,^,  are 
evaluated  using  the  material  properties  at  the  grid  point  at  which 
the  flux  is  being  calculated.  The  data  necessary  for  the  solution 
of  the  flux  at  a  point  is  summarized  in  Tables  X  and  H/ 


Table  11. 

i_ 

JL 

Data,  for  the  S  approximation 

b         c       d 
JL         i       i 

•i 

P 

i 

0 

-1 

-1        0 

— 

1 

♦0.3183 

1 

0 

-1/3      h/3 

-2/3 

2 

*0.361i3 

2 

1 

2/3      V3 

1/3 

2 

♦0.3183 
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Table  IT. 


Data  for  Gauss  quadrative  approximation 


f_ 
0 

1 
2 
3 


-0.8611363 
-0.3399801 
+0.3399801 
+0.8611363 


-1.691*1*625 
-0.319726ii 
♦0.319728U 
+1.691*1*625 


♦1.96770^ 
♦1.06331*0 
♦1.0633110 
+1.96770^ 


+0.31*7851*8 

♦0.6^211*52 
♦0.65211*52 
+0.31*7851*8 


Not  all  pairs  of  values  for  i  and  f  are  relevant  for  each  of  the 

cases  A  and  B  as  a  consequence  of  the  defining  inequalities  for  these 

cases.  The  pairs  of  relevant  values  are  discussed  in  connection  with 

Figure  II.  A  general  point  may  be  defined  to  be  any  grid  point  of  the 

reactor  other  than  a  point  at  the  center,  on  the  centerline  or  in  the 

midplane  of  the  reactor.  A  point  at  these  other  places  is  called  an 

origin  point,  a  centeri»line  point  or  a  midplane  point,  respectively. 

For  Case  A  for  which  |b.|  <|h  land  in  the  case  of  the  first  quadrant 

flux-  there  is  a  further  requirement  (Table  I.  Appendix  A)  that  b 

1 

and  H  have  negative  signs.  In  this  quadrant  there  are  three  Case  A 
fjj.ux  equations,  one  for  each  of  the  pairs  (i«0,  f-0),  (i«l,f«0),  and 
(i-l,f«l).  Similarly,  for  Case  B,  the  requirement  b  I  >|h  leadf  to 
only  one  relevant  combination. 

The  applicable  combinations  and  those  for  the  remaining  three 
quadrants  are  summarized  in  Table  HZ. 
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Fignre  II 

Flux  Configuration 
for  *  particular  energy  group 
at  a  general  grid  point 


m+1 


ra~l 
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Table  III 


Allowable  combinations  of  i  and  f 


Quadrant 

Case  A 

Case  6 

st 

1  / 

i  -  0,  f  •  0 
i  -  1,  f  «  0 
i  «  1,  f  «  1 

i  «  0,  f  =  1 

2nd 

i  -  2,  f  -  0 

i  -  2,  f  -  1 

3rd 

I  -  2,  t  ■  | 

i  -  2,  f  «2 

tf 

i  »  0,  f  -  3 

i  -  1,  f  -  3 
i  =  1,  f  «  2 

i  -  0,  f  -  2 

The  configuration  of  the  flux  at  a  grid  point  which  results  from  the 
allowable  combinations  of  i  and  f  are  shown  in  Figure  II*  If  the  flux 
is  assumed  to  be1  distributed  in^+1  energy  groups  there  will  be  12  (Jc  4tL) 
flux  equation  associated  with  every  grid  point*  For  the  case  considered., 
there  are  three  energy  groups  and,  therefore,  thirty-six  flux  equations 

for  each  point* 

[7] 

The  following  boundary  conditions  are  assumed. 

1*  Continuity  of  flux  exists  at  the  center  line  and  at  the 

midplane  o£  the  reactor,  and  further  the  flux  is  symmetrical 
about  these  boundaries*  This  permits  the  calculation  of 
the  flux  in  only  one  quarter  of  the  reactor*  Figure  III 
shows  the  flux  distribution  at  these  various  special  points* 
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Figure  III 

S  approximation  to  the  vector  flux  distribution 
at  the  boundary  points  in  a  reactor 


HidpJane 

1  '  i\ 


2e  A  second  boundary  condition  is  that  at  the  outside 
boundaries  of  the  reactor  there  exists  no  inward  flux* 
In  order  to  satisfy  this  requirements,  an  extra  set  of  grid  points  is 
placed  Just  outside  these  boundaries ,  and  at  these  points  the  vector 
flux  whose  argument  describes  inwardly  moving  neutrons  is  assumed  to  be 
zero.  Note,  however,  that  leakage  does  exist  at  these  fictitious 
points  and  at  the  outside  boundary  points  of  the  reactor. 
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Nomenclature 


Values  obtained  by  differencing: 'jmd;  integration:  of 
the  transport  equation  overO^  9   and  defined  by 
Equation  (74  to  7d)  of  Appendix  A 

Values  depending  upon  the  quantized  values  of  |i  and 
defined  by  Equation  (7e  and  7f)of  Appendix  A, 

eL  Values  obtained  by  the  use  of  the  Gauss  quadrature 

approximation  to  sum  the  flux  over  n.jpaq»%;e 
defined  by  Equation  (13)  Appendix  A„ 

P  Values  obtained  by  the  use  of  the  S  approximation 

to  sum  the  flux  over  f)   space,  defined  by  Equation  (12) 
Appendix  A* 


r  mean  radius  between  two  adjacent  grid  points 

h 


Neutron  source 

Fraction  of  fission  neutrons  which  are  produced  in  a 
particular  energy  group  -*• 

Number  of  neutrons  produced  per  fission* 


3U 


Ik  , 

T  The  probability  per  unit  length  of  neutron  travel 

that  a  neutron  is  scattered  in  group  k  and  that  a  neutron 

emerges  in  group  A   from  the  scattering  event* 

2  Macroscopic  fission  cross  section 

2  Total  Macroscopic  cross  section 

0(r,z,ii,>Vv)     neutron  flux  which  is  classified  according  to  speed* 

v*  colatitude  and  azimuthal  direction*  cosines  |i  and  Y\» 
and  position  r  and  z* 


A  spacing  of  grid  in  the  radial  direction* 


A  spacing  of  grid  in  the  axial  direction* 

m 


|i  cosine  of  the  angle  0  between  the  axis  and  the 

direction  in  which  neutrons  move*  See  Figure  I 
Appendix  A* 

fj  Cosine  of  the  angle  \D  •  See  RLgure  I  Appendix  A. 

UJ  See  defining  equation*  page  A-12  Appendix  A* 

LU  See  defining  equation,  page  A-17  Appendix  A. 
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Nomenclature 

Superscripts  and  subscripts 

i  index  labeling  the  quantized  values  of  ^?  • 

f  index  labeling  the  quantized  values  Of  \u 

j  grid  position  index  in  the  radial  direction. 

m  grid  position  index  in  the  axial  direction. 

&  neutron  speed  index* 

r 

k  neutron  speed  index  of  group  in  which  neutron  is 

scattered  • 

n  axial  grid  position  diagonally  opposite  grid  point 

at  which  flux  is  being  calculated* 

h  radial  grid  position* 
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Computer  Introduction 

A  computer  program  is  a  sequence  of  logical  instructions  telling 
the  computer  what  mathematical  operations  are  to  be  performed  on  the 
data.  Therefore,  there  are  two  groups  of  inputs,  the  data  and  the 
instructions.  The  data  fed  into  the  computer  are  stored  in  the  fast, 
core  memory  if  there  is  sufficient  room.  When  the  space  is  not 
adequate,  some  of  the  data  are  stored  in  magnetic  drums.  The  only 
difference  in  the  two  storage  processes  is  that  the  data  in  the 
drum  takes  longer  to  locate  and  bring  into  the  core  memory.  Each 
number  or  piece  of  data  is  stored  in  a  register  or  cell  the  location  of 
which  is  designated  by  a  number  called  the  address  of  that  register. 
The  instructions  are  stored  in  sequence  and  are  performed  in  the  order 
in  which  they  are  stored  unless  the  program  directs  otherwise.  No 
unused  registers  are  permitted  to  interrupt  the  storage  as  these  will 
stop  the  computer. 

The  instructions  given  to  the  computer  consist  of  two  parts,  the 
operation  and  the  address  sections.  The  operation  section  is  first 
and  uses  letters  to  tell  the  computer  the  mathematical  operation  to 
be  performed.  The  address  section  gives  the  address  of  the  register 
containing  the  data  on  which  the  operation  is  to  be  performed.  The 
instructions  recognized  by  Whirlwind  I  are  listed  in  the  glossary. 

For  Whirlwind  I  there  are  many  useful  devices  that  make  the 
programming  much  simpler.  A  few  of  these  devices  used  in  this 
program  are  floating  addresses,  preset  parameters,  and  subroutines. 

Probably  the  most  useful  device  is  the  floating  address.  This 
floating  address  or  "flad"  is  used  to  eliminate  the  necessity  of 
referring  to  the  absolute  register  address  in  the  core  memory  when 

37 


writing  the  program.  The  absolute  register  address  of  any  part  of 
the  program  must  be  known  primarily  so  that  the  initial  or  final 
instruction  of  data  is  not  placed  in  a  register  containing  part  of 
the  subroutines  fixed  in  the  memory  of  the  computer* 

In  the  program  key  instructions  on  which  mathematical  operations 
are  to  be  performed  or  which  will  be  referred  to  later  in  the  program 
are  given  floating  addresses*  These  addresses  consist  of  a  letter  and 
a  number s   for  example  al,  a2,  etc.  Any  letter  except  o  and  ell 
may  be  used  and  the  sum  of  all  the  numbers  following  these  letters  may 
not  exceed  23>f>e  When  the  program  is  read  into  the  computer,  it  will 
assign  absolute  addresses  or  register  numbers  to  these  fladSo  A 
table  of  the  floating  addresses  is  produced  by  the  computer  with  the 
solution  unless  it  is  specifically  suppressed  in  the  computer*  ^his 
table  gives  the  absolute  address  of  every  assigned  floating  address 
and  is  particularly  helpful  in  Rouble  shooting  a  faulty  program. 

Another  useful  device  is  the  preset  parameter.  These  are 
numbers  other  than  data  that  are  used  in  the  program.  Preset 
parameters  are  frequently  used  to  designate  the  number  of  times  that 
a  calculation  is  to  be  performed.  Preset  parameters  must  be  assigned 
a  value  before  they  are  used5  and  they  will  retain  this  value  until  it 
is  specifically  changed.  Preset  parameters  are  designated  by  two  letters 
and  a  number,  ^he  first  letter  must  be  a  pj>  u  or  z  and  the  second 
letter  and  the  number  may  be  anything^  except  the  letters  o  and  el  may 
not  be  usedo 

A  third  useful  device  is  the  subrouting.  Its  usefulness  comes 
from  the  fact  that  it  is  a  complete  entity  in  itself  and  may  be  construe^ 
ted  and  trouble  shot  independently.  *or  the  build  up  of  a  complete 
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program  it  is  possible  to  start  with  one  of  the  most  basic  mathe- 
matieal  operational  sequences  that  will  be  performed  repeatedly, 
program  it,  perform  the  hand  calculations  of  the  sequence  and 
then  place  the  routine  on  the  computer.  Modification  of  the 
subroutine  or  correction  of  the  hand  calculations  must  continue 
until  the  results  of  both  are  in  agreemente 

In  arranging  these  subroutines  in  the  main  program  extreme  care 
must  be  used.  The  main  program  must  be  interrupted  to  call  the  sub- 
routine into  use  and  after  the  subroutine  has  performed  its  function 
it  must  direct  the  computer  back  to  the  proper  place  in  the  main  program. 

An  example  of  a  subroutine  is  as  follows I 

Main  program 


isu  p  10 
tl5  isp  x  k 

imr  a  $ 


take  the  next  instruction  from  register 
x  k  and  continue  from  there 

forms  product  of  (°m2)  and  (a£) 


i  STOP 
xks     ita  x  $ 

ica  m  2; 
its  a  7] 


end  of  program 

transfer  into  the  address  section  of  the 
instruction  in  register  x5,  the  address 
that  is  one  more  than  the  address  tl« 

places  m  2  in  a  7 
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ics  m  2 
its  a  8 


places  -m  2  in  a  8 


x59  isp  0  returns  to  main  program  one  register 

after  t  Ij  i.e,;  register  containing 
imr  a  J>* 

Since  the  numbers  used  in  this  program  were  over  four  decimal  digits 
in  lengthy  floating  point  arithmetic  was  required*  In  this  system  two 
registers  are  required  to  express  a  number*  Each  of  these  registers  is 
capable  of  containing  16  bits  of  binary  number  information*  The  first 
bit  of  the  first  register  expresses  the  sign  of  the  number  and  the 
first  bit  of  the  second  register  expresses  the  sign  of  the  exponent  of 
20  Six  bits  are  required  to  locate  the  binary  point  of  the  number 
leaving  2k  bits  to  express  the  actual  number*  In  certain  stages  of  the 
program  it  is  desirable  to  handle  single  length  register  numbers*  The 
instructions  for  the  handling  of  these  numbers  are  similar  to  those  for 
double  length  registers  except  that  the  i  does  not  precede  the  instruc- 
tion* When  operating  on  double  length  registers  the  computer  is  said 
to  be  in  the  interpretive  mode,  when  on  single  length  registers  it  is  in 
the  Whirlwind  mode*  To  enter  the  interpretive  mode  from  the  Whirlwind 
mode  it  is  necessary  to  use  the  instruction  "IN"  and  to  reverse  the 
direction  it  is  necessary  to  use  the  instruction  "OUT"* 
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GLOSSARY  OF  INSTRUCTIONS 

TERM  DEFINITION 

ica  x  clear  MRA  and  add  contents  of 

register  x. 

ics  x  clear  MRA  and  subtract  contents 

of  register  x. 

iad  x  add  contents  of  register  x  to 

what  is  in  MRA. 

isu  x  subtract  contents  of  register  x 

from  what  is  in  MRA. 

imr  x  multiply  what  is  in  MRA  by  the 

contents  of  register  x  and  round 
off  product  to  fifteen  digits* 

idv  x  divide  contents  of  MRA  by  the 

contents  of  register  x. 

its  x  transfer  the  contents  of  the  MRA 

to  register  x  losing  what  was 
previously  in  register  x.  The 
MRA  remains  unchanged. 

iex  x  exchange  the  contents  of  the  MRA 

with  the  contents  of  register  x« 

isp  x  take  the  next  instruction  from 

register  x.  Does  not  affect  the 
contents  of  the  MRA:*  An  uncon~ 
ditional  instruction  used  to  break 
the  sequence  of  operations* 
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Glossary  of  Instructions  (Continued) 

TERM  DEFINITION 

icp  x  if  the  contents  of  the  MRA  are 

negative^  take  the  next  instruc- 
tion from  register  x  and  continue 
from  there o  If  the  contents  of 
the  MRA  are  positive,  take  the 
next  instruction  in  sequence*  A 
conditional  instruction  used  to 
break  the  sequence  of  operations. 

isc  x  select  counter  number  x.  Without 

this  instruction  counter  zero 
will  be  used  whenever  +C  appears. 

icr  x  cycle  reseto  Sets  index  register 

of  counter  to  zero  and  criterion 
register  to  x© 

ict  x  cycle  transfer  Increases  con- 

tents of  index  register  by  one* 
If  contents  of  index  register  greater 
or  equal  to  contents  of  criterion 
register!,  set  index  register  to 
zero  and  do  next  instruction 
in  sequence o  If  contents  of 
index  register  less  than  contents 
of  criterion  register^  take  next 
instruction  from  register  x« 

ici  x  increase  the  contents  of  index 

register  by  number  x» 

icd  x  decrease  the  contents  of  index 

register  by  number  x. 
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Glossary  of  Instructions  (Continued) 

TERM  DEFINITION 

iti  x  transfer  the  right  eleven  digits 

of  the  index  register  into  the 
right  eleven  digits  of  register  x© 

iat  x  add  contents  of  index  register  to 

the  contents  of  register  x  and 
store  the  result  in  the  index 
register  and  register  x© 

ita  x  replace  the  address  section  of  the 

instruction  in  register  x  with  the 
address  that  is  one  more  than  the 
address  of  the  register  containing 
the  last  isp  or  icp  if  the  contents 
of  the  MRA  are  negative© 

iTOA  record  the  contents  of  MRA  on 

direct  printer  (typewriter)* 

iMOA  record  the  contents  of  MRA  on 

magnetic  tape  for  delayed  printing. 

iSOA  record  the  contents  of  MRA  on 

oscilloscope  for  photographing© 

iFOR  this  instruction  provides  an 

automatic  device  for  obtaining  a 
suitable  layout  of  output  data 
in  columns s   lines  or  blocks© 

Notes  The  ibove  output  instructions  are  usually  followed  by  a  series 
of  letters  and  numbers  that  indicate  the  desired  form  and 
arrangement  of  the  output© 
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Computer  Logic 

The  program  starts  at  1  in  the  main  program  by  entering  the 
data,  program  and  preset  parameters.  The  initial  calculations  are 
•with  the  reactor  radius  equal  to  the  -value  placed  in  register  rli. 
With  this  radius  the  program  then  shifts  to  the  flux  iteration  and 
convergence  test,  ho     In  this  routine  parameters  are  set  to  control 
the  number  of  iterations  allowed  at  one  radius  and  to  traverse  the 
computer  thru  all  points  in  the  reactor.  With  a  value  of  j  and  m, 
defining  a  point  in  the  reactor,  the  program  shifts  to  Flux,  li.J>, 
where  the  addresses  on  the  drum  of  the  required  data  are  calculated, 
the  data  entered  in  the  fast  memory  and  the  routine  appropriate  to 
the  point  is  selected  and  control  is  shifted  to  that  routine© 

The  first  calculation  at  any  point  is  the  scalar  flujie  for  all 
three  groups  and  the  second  is  the  input  from  fission  and  scattering 
from  an  upper  energy  group*  The  first  point  to  be  calculated  in  each 
cycle  is  the  origin0  This  is  the  point  of  normalization  of  the 
inputs  to  the  fast  groupo  The  normalization  constant50  9   is  defined 
as  the  reciprocal  of  the  fast  group  input  at  the  origin.  With  this 
constant  all  other  inputs  (all  groups  and  all  points)  are  normalized* 

As  explained  in  the  previous  section  of  the  procedure  certain 
quadrants  of  the  flux  are  required  for  each  point  depending  on  its 
location©  Prior  to  sequencing  to  another  point  all  fluxes  appropriate 
are  calculated* 
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T$b  fast  source  input  is  then  compared  to  the  previous  fast 
source  input.  If  the  new  value  is  not  within  a  set  specified  fractional 
value  of  the  previous  value  a  negative  quantity  is  introduced  into  s  10* 

When  the  entire  grid  has  traversed  a  check  is  made  to  see  if  the 
maximum  allowable  attempts  at  source  convergence  have  been  made*  If  this 
maximum  number  has  been  made  the  computer  will  stop*  If  not  a  check 
will  be  made  of  switch  s  10^  As  indicated  in  routine  b«  another  •vale 
will  be  made  thru  the  reactor  if  a  10  is  negative,  and  if  s   10  is  positive 
control  will  return  to  the  main  program  to  repeat  the  above  procedure 
based  on  the  radius  placed  in  r  3« 

When  a  gamma  is  calculated  for  this  radius^  control  is  again 
returned  to  the  main  program  but  at  a  different  pointo  All  subsequent 
returns  are  at  this  point*  A  check  is  made  to  determine  if  the 
maximum  allowable  number  of  iterations  of  radii  have  been  attempted© 

A  test  for  criticality  is  made  to  see  if  gamma  is  within  a 
specified  tolerance  of  unity*  If  it  is  note  a  new  radius  is  found  by 
extrapolation  or  interpolation  based  on  the  most  recent  values  of  R  and 

Xand  the  best  previous  value  of  R  and  ft  o     In  this  determination  of 
radius. the  gamma  that  is  closer  to  unity  and  the  associated  radius  are 
preserved  for  use  in  later  extrapolation*  The  previously  outlined  cal- 
culations are  repeated  with  this  radius*  If  o  lies  within  the  specified 
tolerance  of  unity  control  remains  in  the  main  routine  and  the  print-out 
routine*  which  has  been  stored  in  the  drum*  is  entered  into  the  fast 
memory* 
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The  computer  is  then  sequenced  thru  the  grid  calculating  and 
printing  out  scalar  fluxes  at  each  point  for  each  group* 

If  the  flux  pattern  is  not  as  desired  the  properties  of  the  regions 
are  changed  and  a  new  tape  of  data  cut.  The  program  is  then  run  again  to 
determine  the  flux  pattern* 
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RESULTS 

Numerous  runs  of  subroutines  and  of  the  entire  program  were  made 
on  Whirlwind  with  the  program  contained  in  Appendix  B.  This  program 
was  adapted  for  a  reactor  radius  of  75  and  95  centimeters  and  a  height 
of  168  centimeters •  The  grid  placed  on  the  reactor  was  a  6  by  6  grid 
giving  for  the  first  trial  radius  a  grid  size  15  cm  by  31*  cnu 

The  final  run  had  programmed  a  first  trial  radius  of  75  cmj  the 
maximum  allowable  number  of  iterations  with  a  constant  radius  was 
twenty  and  the  tolerance  of  convergence  of  the  fast  source  input  was 
three  percent.  The  computer  stopped,  as  programmed,  upon  completion 
of  twenty  iterations  at  a  set  radius  when  the  fast  neutron  input  at 
all  points  was  not  with  the  tolerance  specified*  The  computer  running 
time  per  iteration  is  two  minutes.  The  computed  values  of  scalar 
flux  at  each  point  and  energy  g^oup  were  printed  out  for  each 
iteration*,  The  normalizing  constant  decreased  from  0#27  for  the 
first  iteration  to  0o01  for  the  twentieth*  Some  values  of  the  scalar 
flux  at  the  following  point  were  negative 

a)  interface  of  core  and  reflector « 

b)  reflector 

c)  outside  boundary  of  reflector* 
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DISCUSSION  OF  RESULTS 

The  negative  values  computed  for  scalar  fluxes  are  known  to  be 
unrealistic.  The  possible  causes  of  this  error  are 

1)  improper  programming  of  the  theoretically  derived 
flux  equation, 

2)  grid  size  being  excessive* 

The  theoretical  background  for  the  flux  equations  is  sound,  and 

these  equations  have  been  programmed  for  other  geometrical  arrangements 

[11]* 

on  other  computers  and  give  correct  flux  distributions    •  The 

need  for  a  complete  hand  calculation  is  well  realized,  but  due  to  lack 
of  time  these  hand  calculations  were  not  completed  fully.  Therefore, 
there  is  no  reason  to  assume  or  expect  that  these  equations  have  been 
programmed  correctly.  To  accomplish  the  required  hand  calculations, 
computer  verification  of  the  hand  calculations,  and  program  trouble 
shooting  it  is  estimated  that  another  three  to  four  weeks  full  time 
would  be  required© 

For  accuracy  in  evaluation^ of  neutron  transport  equations,  the 
grid  spacing  should  be  less  than  or  approximately  equal  to  the  mean 
free  path  of  the  neutron*  If  this  condition  does  not  prevail,  the 
accuracy  and  reliability  of  the  calculations  decrease  rapidly. 

A  two  fold  increase  in  the  number  of  grid  points  will  increase  the 
time  of  computation  by  a  factor  of  eight.  Any  increase  in  the  time 
required  for  one  iteration  thru  the  reactor  is  not  allowable  since  ndir 
two  minutes  are  required. 
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The  code  developed  is  in  the  •interpretative'  or  'floating 
point'  mode  for  "Whirlwind,  for  the  reason  that  this  mode  is  simple 
to  learn  and  to  trouble  shoot*  Whirlwind  has  also  available  a 
•Whirlwind1  or  'fixed  point'  mode  in  which  the  speed  is  increased 
by  a  factor  of  roughly  3S>»  Use  of  the  fixed  point  mode  would 
make  the  two  dimensional  S  calculation  on  the  Whirlwind  computer 
practical* 


.  60 


CONCHJSIONS 

1)  The  code  as  presently  programmed  is  incorrect* 

2)  The  grid  spacing  as  used  is  much  too  coarse*  This 
difficulty  is  not  inherent  in  the  coded  program* 

3)  The  mode  of  computer  operation  is  unsuitable  for  an 
iterative  process  such  as  this  because  it  results  in 
a  code  that  progresses  too  slowly* 
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RECOMMENDATIONS 

1)  Prior  to  the  use  of  any  computer  a  formal  or  semi^tforaal 
course  in  programming  should  be  completed* 

2)  A  grid  spacing  comparable  to  the  mean  free  path  of  the 
neutron  should  be  selected* 

3)  The  fastest  possible  mode  of  calculation  on  the  computer 
available  should  be  used* 

k)   For  a  coded  program  of  this  difficulty  a  faster  machine 
with  a  much  larger  fast  memory  is  very  desirable* 
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Appendix  A 

Development  of  Equations  Using  Sp  Approximation  to  the 
Boltzman  Equation  for  Cylindrical  Geometry 
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Equations  will  be  developed  for  the  time  independent  case  using 
cylindrical  geometry. 

The  vector  flux  at  a  point  in  space  is  a  function  of  two  variables^ 
one  describes  distance  along  the  vertical  axis  of  the  cylinder,  and  the 
other  describes  radial  distances.  The  axial  variable  will  be  designated 
as  |x  and  the  radial  variable  as  n  •  From  Figure  I  it  can  be  seen  thatt 
n-  «  cos  V 


p.  =  cos  0  and 

sin©"  \1  -  [i 

sin  y  «  \|l-^2 

and  where SX   represents  the  unit 

vector  in  which  direction  neutrons 

travel  at  a  pointo  Therefore , 

V  0   «  i  — ^  +  i  —  -rfr;  *   i  — c 
r  i>r    Qr  3Y    z  <>z 


Figure  I 

Components  of  neutron  flux 
at  a  space  point. 


and  -A-  *  i  sin  0  cos  1 
r 

i  d  0  ♦  a 

9  sin  9  sinV  "rr   i*  c©a  0, 


Thersfora, 


A 
-A.  e 


s  sin  9  cos 


sin  9  sin 


*   * 


2i 


NOb  ♦  cos  9 
*r  6% 


If  we  substitute  the  expression  for  yl  and  \i  which  were  derived  from 
Figure  I  in  equation  lj>  we  find 


jz  •  Vi 


"21  ii  d-n  jjt""?  >0 


(i) 


3  z 

The  quantity  J£ e  0  represents  the  leakage  of  neutrons  out  of  an 
element  of  wlume  per  unit  time  in  the  direction  xk   P«r  unit  solid  angle 


(2) 


A=l 


In  the  %±me   independent  case  this  leakage  pltra  the  rate  of  loss  of  neutrons. 


/E  09   per  unit  time  must  equal  the  source  strength^  S„  In  problems 

in  which  cross  sections  vary  with  energy  and  neutrons  alter  their  speeds, 

[6] 

the  multigroup  approximation  of  the  Boltzmann  equation  is  introduced* 


Vf. 


d 


2 


J± 


*ll± 


T 


i  *75?  az"1I7   "fi^1 


i. 


where   2  is  the  total  cross  section  within  the  energy  group  UC   and 

S  is  the  neutron  source  term  which  results  from  fission  and  energy 

degradation  of  the  total  vector  flux  at.  the  point.  Slowing  down  is 

incorporated  into  the  source  term. 

Equation  3  is  identical  in  form  to  the  time  dependent  multigroup 

equation  for  spherical  geometry  in  which  the  z  term  is  replaced  by  a 

term  representing  the  rate  of  change  of  flux  with  respect  to  time. 

The  S  approximation  is  a  method  of  evaluating  integrals  over  a 
n 

variable  in  which  the  variable  is  quantized  into  n  intervals  and 
intermediate  values  of  a  function  of  the  variable  are  approximated  by 
a  straight  line  between  points  at  the  boundaries  of  the  intervalso 
We  shall  consider  v?  space  as  the  first  variable  which  will  be  Inte- 
gra tedo  The  flux  at  any  point  Figure  II 


(r5Zc»n  )  in  Figure  II  can  be 
represented  by  the  following 
equation  where  i  represents 
the  interval  number© 


9(^^»w    Actual 
flux 


7(    Mesh 


approximation 


(3) 
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<*s*s> 


,      Alli-i,,  ,      VI       rff 


(« 


Foz»  the  S     approximation  i  ■  +lfl  05  and  ~lo     If  we  denote  10    -  n  . 
as  Ao   and  0(rpzam»o  .  )  as  JT  we  obtain 


If  we  substitute  this  in  equation  3*  we  find 


A  ^  r  Ao  "^  r  r  A.  <5>n 

i  i  i  L 

^ (1~  i  )(ii-i}  3V°X  >  W    ?V  ^  ^r^   ^V4 


C3  S 


^VU^  ^ 


li"V 


A„    ^  1-p.' 


Since  the  source^  g,,  is  isotropic  integrating  this  expression  over 
from  f)  ton    s  we  obtain  n  equations  in  the  form 


OI 


i     c)  r  i    r  "^  z 

a* 


3  V"      La         3 

*  h    ^    ^     ♦  gZ  £T  +  d 


i       £  r 


r 


«±  g  s 


where 


!ti*1i-i 


i  v 


(7«) 
(7b) 
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4|  !*•  z% 


i      3  y  i  .  |?  ' 

•i"2  (W)     gj^-Ljjp,.  (7f) 

A  different  method  of  finite  differencing  is  used  for  the  other 
directional  variable5*4f9  since  use  of  an  expression  analogous  to  (1;) 
for  n  in  the  integration  of  the  equation  (5)  with  respect  to  |i  would 
result  in  an  extremely  complex  expression.  Instead,  the  equation  (5) 
is  directly  quantized  with  respect  to  \i  giving  the  result^ 

\hT~  +\     r     +hf"tT"  +8f2^    ' 

(8) 

i<Jr  i        r  f     $  z  f  if 

We  may  now  proceed  to  evaluate  the  source  integral  at  a  point© 
For  isotropic  scattering  and  three  energy  groups  of  neutrons 


S(r^)  -  yv    Jx  2f  (l(r9*)  *  J^  VV(^z). 


(9) 


where 


k      x     ■'*        +1         dfl 


As  in  equation  (ii)  f(vD%9ntA)   is  approximated  by 


'  A    Mr^u)  ♦  :     PT  (r^ii), 


^#»*i*»  )  s 


A-fc 


and 


>5z3|i) 


2^  JL- 


Z 


^  A  ^  4  «     I 


+  0 


i-1 


I1 


^  -       -7 


(10) 


Upon  integration 
0(r«,z,|i)  = 


2.  r  4      | 

-S5  i~FC  -vnT+Visin"\i-i 
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i.isin"\3] 


i 


n  .     sin     /j 


2 


<  10         sin     n 


(11) 


u       i 


i-1 


n 


««n..   m-r^    -U-^ 


If  values  of  fl  o  *  ♦!<,  0  and  =1  (S  approximation)  are  substituted  in 
the  above  expression,  we  find  that 


fyr,i,|0  -  ^   Jo  P. 


(r^ji) 


(12) 


where  P  -  0„3183*  P  -  0o363li<>  and  Po  -  0o3l83  . 
o  1  2 

[g] 

By  usa  of  the  Gauss  quadraiive  method   the  remaining  integral 

can  be  evaluated  in  the  form  of  a  sum  of  the  products  of  coefficient 
times  the  vector  flux  where  p.    values  are  the  n+2  positive  roots  of 


A-$ 


The  table  below  lists  the  values  of  u  and  the  coefficient  H  used 


f 

Mf 

H 
f 

0 

-0.8611363 

0.3U78S1|8 

1 

-0. 3399810 

0o6^211i52 

2 

+0* 3399810 

Oo  6£212£2 

3 

♦Oo 8611363 

0.31*7851*8 

Therefore g  since  H  =  H  and  H  =  H  the  scalar  flux  at  a  point  is 
5       o    3     1    2  * 

k„,   %   1  _    2  _  ki,r 


Z)  "floHf  i2oPi  V'r(r,.), 


(13) 


and  the  source  term  equals 


-jri  ,  1  \    V  (r,a)  ♦  ,  I      X\  fyr.z),  .  (11,) 


5  k«l  ^f   "  vr*z;  *  k=l 


where 

p  is  the  propability  that  a  fission  neutron  will  have  an  energy jt  and 

Jtk 

T^is  the  probability  that  a  neutron  having  an  energy  k  will  be 

scattered  into  the  energy  group  X  5  k>  £* 

The  scattering  probability  requires  further  explanation..  We  define 

th 

\l  g      as  the  probability  that  a  neutron  which  is  scattered  in  the  k 
/.k 

„tfa 

group  will  land  in  the  it       group^  k  as  the  energy  group  in  which  th© 

scattering  occurs  (l  for  fastr,  2  for  epi thermal  and  3  for  thermal)  and 
Si   as  the  energy  group  into  which  the  neutron  is  scatteredo  The 

energy  range  covered  by  scattered  neutrons  is  E»aEp  and  the  proba- 

AE 

bility  of  landing  in  the  interval  AE  is  _  '  rri  <>     Let  E  be  the  energy 

E— ctE 

of  the  neutron  before  scattering© 

E   <  E  <  E  o 
k+1      k 

A-6 


The  neutron  flux  per  unit  energy  interval  is 

0  E 

d  0  ' 

Then  ,=~£'  dE  £  (E)  is  the  number  of  neutrons  scattered  per  unit  volume 

"       8  [9j 

and  per  unit  time,  but 


dE   2  E  | 

d#  q  dE 

Therefore,,  -±  dE  S  (E)  -  7  ■—  (l£) 

dE      s      ?   E 

E„  -  E  . 

— —7" — is  the  probability  a  neutron  will  land  in  the  /  th  group* 

E(l=>cO 

y  is  the  number  of  neutrons  scattered  in  the  k  th  group  that  land 
in  the  X  th  group©  Since  the  absorption,  sources  and  leakage  in  a 
large  thermal  reactor  are  small,  the  slowing  down  density  qs  is  rather 
constant  over  an  energy  group<>  This  fact  is  used  to  compute  the  num= 
ber  of  neutrons  scattered  from  the  group  k  to  groupyJ^  • 
The  equations  for  each  case  are  listed  belowo 

Case  Equation 

lo  E,  >  aE  y  B     x  r     ,      Z.       n  (E  -E   )(■=*—  -  -*-)  (16a) 

£+1        k  T       (l~a)  [ln(E  /E       )  K   .  *1  AE  ^       * 

k  k*l  k+1   k 


2.  E  -  >  a  E.  jtk  s 


'I  *  u  k 


aE  >E,   >aE 
k  jt+1   k+1 


r = 


(1-a)  ln(E  /E  ,) 
k'  k+1 


A=7 


E      "  Q  "  E 
.  kl1         k 


(16b) 


3.  l^>«Bk 

J?+l  k+1 


E         "  E 
k+1        k 


X  a    In 


k+1 


(16c) 


t a  \+1  <  V  <  a\ 


E.      <  a  E 
X  +1  k+1 


Zk 


(I6d) 


r   '  (1-a)     In 


Sv3 


E 


E 


E 


_  k+1        Ji 


-r-aln 
E 


a  E 


k+1 


£„  E  .  <  a  E 

JL  k+1 


Zk 


r 


(I6e) 


For  the  cases  where  the  scattering  occurs  within  the  energy  group 
we  find 


E  .      ,  >  a  E  . 
JL+%  A 


E.       <  a  E 
Ml  JL 


JLt 


r  = 


l~a 


r       e       -  e 

L   *ejl  lnE^+l 


(I6f) 


JLfi 


L- 


x 


In     E,  /E 

Z  '  ^  +1 


We  may  now  proceed  to  evaluate  equation  k  by  integrating  over  the 
r  and  z  variables <>  However,,  before  doing  so  we  must  consider  the  direction 
in  which  integration  should  be  carried  out.  It  is  always  best  to  integrate 
in  the  direction  of  neutron  flow  since  any  errors  produced  in  the  process 
of  integration  will  be  diminished  as  we  proceed  further  in  this  direction* 


A-8 


If  integration  is  carried  out  in  a  direction  opposite  to  that  of  neutron 

[7] 

flow,,  any  errors  produced  will  be  magnified. 

We  must  now  difference  our  equations  with  respect  to  r  and  z<> 
We  can  arrange  terms  in  equation  (8)  in  order  to  facilitate  the  integra- 
tion0  Thus  equation  (8)  becomes 


f 


g   ^s 

if 


(17) 


We  must  now  further  divide  the  results  of  the  integration  into  two 
separate  cases* 

Case  A  ■=>  The  component  of  velocity  in  the  r  direction  is  less  than 

the  component  of  velocity  in  the  z  direction* 

Case  B  -  The  component  of  velocity  in  the  r  direction  is  greater 

than  the  component  of  velocity  in  the  z  direction© 
If  we  consider  Case  A5  and  difference  r  and  z  as  illustrated  in  the  figure 
below5  we  obtain  four  distinct  subcases©  These  will  be  classified  in 

accordance  with  the  direction  of 
travel  being  determined  by  the  signs 


Figure  HI 
Position  Mesh 


M+l 


of  b     and  h  •     A  positive  sign 
i  f 

indicates  a  direction  of  travel  away 
from  the  center  of  the  core*  The 
various  cases  are  listed  in  Table  I» 
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Tabl«  I 

In  all  A  cases  b  |<|h  I 
1  i'  '  f  I 

In  all  B  cases | b  |  >|h  | 


Subcase 


Ml,   Bl 


A2,  B2 


A3a  B3 
M9  Bl* 


+ 


Considering  subcase  A-3  and  integrating  equation  £  and  knowing 


that 


r.  2  L  j5m      .i-l 


0-l*n 


and 


/  pCr^z)  da  -  — 

m-1 


r 


/A,f 


l«f 


j5m     j5» 


where  A  =  r  -  r 

j    j    j-1  and  A  «  z  -  z 

m    ra    m=l 

we  obtain  the  following  equation 
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j-l,m-l 


ei  1 
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j-1 
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j^m  j-l,m-l 

2 
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g. 


g       V1 
gf  j-1 

^-i  _  JL         Jla  :KUm-l 


J-l       r 


j-l 


b    -  4 
i      vi 


^-±'r  *  ,V1,r .  V"1,r  -  "V"1** 


j,m  J5m~l  J~l*m 


j*»l,m-l 


i  f 
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S         +      ^S 


j-l,m-lj 


(18) 


Where  Z   is  the  total  cross  section  of  the  material  at  the  center  of 

3-1 
the  grid  in  figure  k$   i.e«a  the  material  to  the  left  and  below  the  ,j5m 

—      r  ♦  r 
point  and  V        is  j    jKL 

d"1      2 

The  radial  direction  of  integration  is  determined  by  the  sign  of  fj  9 

ioeo5  by  the  sign  of  the  b  or  d  •  Consequently,  the  signs  of  Y)   9   b  , 

daS>   and  A  all  change  together.  Likewise,  the  direction  of  integration 
<j 

axially  is  determined  by  the  sign  of  u.5  i.e.,  by  the  sign  of  h  o 

Consequently  the  signs  of  p.,  h  and  A  all  change  together.  If  we 

f     m 

use  absolute  magnitude  signs  for  these  quantities  all  Case  A  equation* 
will  be  identical  insofar  as  the  signs  of  various  terms  are  concerned© 
Similarly^,  all  Case  B  equations  will  be  identical* 
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If  we  now  add  the  following  expression  to  both  sides  of  equation  (18) 


V**     ♦  V'1*1-   - k  <t'f  -  V"1>f 

j-l,m-l  jfl,ra-l  j,ia  j,m 


we  obtain 

4J>.il 


("v1-  ♦  V"1**  -  v*f  -  V-1'* 

j,m  j-l,m-l  j,m-l  j>m-l_ 


J^jn-1 


'f        j-1      r 


o-i- 


j.irf  j-l,m-l 
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+    r 
r 


j-l 


l2^ 


Z  m-1 

*f       j-1 


j,m  j,m-l 


(18a) 


V1,f  -  V"1,f 

3-l,ja  j-l,ra-l 


ei8f 


V    ♦  V"1,f  -  V,f  -  V_1*f 

j-l,m-l  j-l,m-l  j,m  j,m 


If  this  expression  is  multiplied  by  A 


*U    IaI  ^  ^ 

60     where  A.    co     -  - 


m 


we  obtain 
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a-JS 


V,f  ♦  V_1,f-V,f  -V_1»* 

0,m  j,m  j»m-l  J»i»-1 


OL) 


y      blJJ/g     V"1^ 


-x  •,  M 


k^ 


1  8f  Fj-i 


^f-f    +  A  A)       . 


(18b) 
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Where 


.i1    2    V^r     3-1     § 


Fj-l,m-l 


jj   m-1         °i 


i1       2       \Jf       j-1 


>1/J 


|^l      2     ^f    J-l      r~^/  ">1,»-1        2    '    j\ji,»      j-1,-1 

—  P^"1,f  ♦  V"1'*  -  V"i,f  -  V-1*'  1 

2        /^ j,m  j,m-l  j-l>w  j-l,m-l] 


(18c) 


And  finally  transposing  terms  in  equation  (l8b)  we  obtain  an  expression 
for  the  vector  flux  at  the  point  ;},»• 

i      b  a-Ji'f)  f*V»f  ♦V-1*f-  VMW'V 


,,.yytotA 


OJ   2  I  V>1   r 


j-l 


(lgd) 


By  use  of  similar  procedures  for  developing  the  equations  for  tha 
three  remaining  A  subcases  we  obtain  the  three  equations  on  the 
following  pages* 
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For  subcase  A.-U 
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where 


>k  -Lw 


-kl^fe 


c 

1  i 


j+1        r 


V-f 


d*l 


!i\ 


u.QI  >i  .  A- 


2  Jl\       J,m        $»1,b-M 


FJ+l,n  ^j+l,n-l_J 


V1'* 


j+l,m-l 


(19b) 


A~lit 


For  subcase  A-l 
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For  subcase  A-2 
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Figure  V 

Position  Me$h 
Case  B 


A  similar  procedure  is  used  in  deriving  the  equations  for  the 

8  cases;  except  that  as  shown 
in  Figure  $   the  r  and  z  com- 
ponents of  the  flux  are  different 

because  b  &•  h  •  The  method 
i    f 

of  deriving  the  flux  equation 
for  case  B-3  is  outlined  on 
the  following  page.  The 
equations  for  the  other  B  cases 
are  listed  also. 
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To  both  sides  of  equation  (22a)  we  add 


and  define 
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•if 


I M lb 
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we  obtain 
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Prom  this  expression 
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WAr*  *3 


i  gf 


j-i    o-i 


(22c) 


where  A  \  is  defined  by  equation  (18c) 


For  case  B— 1|. 


0sra  ,      I    |A_.|/    6  ..  ,         c_.     \ 


2*21 I 


b  |JJ 
i      2 


1*  3+1    W 


(23) 


where  A  is  defined  by  equation  (19b) 
k 


For  case  B-l 


3  s* 


Mu-'WV*  ♦  V_1»*-'V'*  ,-  V-1>f  Va 

f     IjTj     "  L      3*1.*  j*l,m        "l+l,m»l        yj*l,m*lJ     1 


2  *f  **   v 


(2U) 


where  A  is  defined  by  equation  (20b) 
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For  case  B-2 


l±9t  j\ 


-1    <°  L   yj-l»n  J-l,»  Vi.m+1    ZjtoJa! 


A2 


li 


\ff  j-a,    5     y  (25) 


where  A  is  defined  by  equation  (21b)« 

The  vector  fluxes  defined  by  equations  (l8d),  (l°a),  (2fca),  (21a), 
(22c),  (23),  (2li),  (25),  are  valid  for  all  values  of  n  with  the  excep- 
tions of  the  end  points  of  our  u-  space  where  \i  ■  +1«  At  these  points 
trouble  is  experienced  with  singularities*  However,  a  modified  equation 
may  be  derived  by  letting  the  terms  containing  r  in  all  equations 
following  equation  (l)  be  replaced  by  zero,  since  for  [i  ■  ♦  1  the; 
term  in\f?  in  equation  (l)  vanishes.  This  special  equation  was  not 
needed  since  values  of  p.  were  chosen  which  avoided  these  singular 
points. 

The  procedures  described  above  result  in  n  equations*  One  more 

equation  is  needed  and  is  provided  by  setting  n  ■  -1  in  equation  (3)» 

The  development  above  may  be  used  when  the  values  of  b,  c,  0- ,  0 

o  mjlt9 

and  £  used  above  are  replaced  by  -1,  0,  0(  *  "-1),  0,  and  1,  respectively* 
These  values  reduce  equation  (6)  to  that  found  for  the  special  case* 
Thus  equation (6c)  becomes  with  i  ■  0 

J,ra  0,f 

2      lgf  0-1 
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where 


+1  - 


k 


Fj-l,m-l       2  I    a 


1.    a. 


s      ♦ 


m-l,m-l 
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Nomenclature 

b^   ]  Values  obtained  by  differencing  and  integration  of  the  transport 
*■   /  equation  over  *  ,  and  defined  by  Equations  7a  to  7d  of  Appendix 

•i  j  *• 

g  ~\       Values  depending  upon  the  quantized  values  of  y,  ,  and  defined 
h      by  Equations  7e  and  7f  of  Appendix  A. 

i      a  unit  vector  in  the  direction  indicated  by  a  subscript. 

9 

^p     Fraction  of  fission  neutrons  which  are  produced  in  a  particular 

energy  group  X  • 

r      The  radial  distance  to  a  point . 

r      The  mean  radius  between  two  adjacent  grid  points. 
h 

z      The  axial  distance. 

E      Neutron  energy  from  a  point*  See  Figure  1,  Appendix  A« 

H      Values  obtained  by  the  use  of  the  Gauss  quadrature  approximation 
to  sum  the  flux  over  p.  space,  defined  by  Equation  13,  Appendix  A* 

P.     Values  obtained  by  the  use  of  the  S  approximation  to  sum  the 
flux  over  |x  space,  defined  by  Equation  12,  Appendix  A* 

oc    Ratio  of  the  minimum  neutron  energy  after  scattering  to  its 
energy  before  scattering* 
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&       Grid  spacing  of  a  mesh* 

Y)  Cosine  of  the  angle  y 

9      Angle  between  the  axis  and  the  direction  in  which  neutrons 
travel.  See  Figure  I,  Appendix  A. 

p,  Cosine  of  the  angle  ©. 

1    th 
H»     The  probability  that  a  neutron  which  is  scattered  in  the  k 

th 
group  will  land  in  the  d   group* 

D     Number  of  neutrons  produced  per  fission. 

2      Total  macroscopic  cross  section* 

2      Fission  macroscopic  cross  section, 
f 

Q.  k 

<v/        The  probability  per  unit  length  of  neutron  fcravel  that  a 

neutron  is  scattered  in  group  k. 

#  (r,z,  Neutron  flux  which  is  classified  according  to  speed,  v, 
r   colatitude  and  azimuthal  direction   cosines  \i  and  a  ,  and 
position  r  and  z* 

T     Colatitude  angle  denoting  neutron  travel  with  respect  to  the 
radius  •  See  Figure  I,  Appendix  A# 

CO    See  defining  equation  Page  A-12,  Appendix  A* 

i,y        See  defining  equation  Page  A-17*  Appendix  A. 

-^   Direction  of  unit  vector  in  which  neutrons  travel* 
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Superscripts  and  Subscripts 

i  Index  labeling  quantized  values  of  rt  • 

f  Index  labeling  quantized  values  of  \l» 

j  grid  position  index  in  the  radial  direction* 

m  grid  position  index  in  the  axial  direction* 

j£  neutron  speed  index, 

k  neutron  speed  index  of  group  in  which  neutron  is  scattered* 
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Appendix  B 
Whirlwind  Program 


87 


Prior  to  the  use  of  this  program  the  following  preset  parameters 

must  be  defined  and/or  calculatedt 

pzl  *  number  of  grid  points  in  a  radial  direction  less  one 

pz2  ■  number  of  grid  points  in  an  axial  direction  less  one 

pz3  ■  number  of  iterations,  N,  allowed  at  a  set  radius,  in 

attempting  to  achieve  a  convergence  of  source  pattern*-  , 

pal;  ■  number  of  adjustments,i.Nf,  of  radiua_Allowed  in 

attempting  to  achieve  criticality. 

pf2  fc  drum  address  of  cross  section 

pfl  +  92(pzl  +  2)  (pzl  +  2) 

pf3  ■  92(pz2  +2) 

pfl*  -  6(pz2  +  1) 

The  preset  parameters  listed  below  were  used  in  the  development  of 

the  program  and  have  the  values  as  indicated 

pfl  «  5°(X)$  drumt  address  of  first  flux  at  origin  point* 

pff?  c  70itj  address  in  fast  memory  at  which  the  print-out 

routine  will  be  placed  when  needed 

pflO  ■  26£0$  address  of  the  print-out  routine  in  the  drum 

pfl5>  ■  3 00 j  length  of  the  print-out  routine 

pfl3  =  litfOj  address  in  the  fast  memory  of  the  vector  flux 

subroutine  of  the  quadrant  being  calculated. 

It  is  necessary  also  to  furnish  the  program  at 

l)  a9,  the  value  of  +(2pzl).0#  This  value  is  used  in  the 

calculation  of  _± 
^l 
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2)  rii,  the  value  of  the  first  trial  radius  of  core 

3)  r3*  the  value  of  the  second  trial  radius  of  core 
For  a  three  group  solution,  values  of  each  of  the  thirty-six 

fluxes  and  the  three  sources  at  each  point  are  assumed  and  placed  on 
tape  in  the  order  indicated  in  Table  I. 

Table  I 
Sequence  of  fluxes  at  a  point 


0 

1 
2 
0 
1 
2 
0 
1 
2 
0 


0  .<  f  <  3     0-  <  i  <  2       0<^<2 
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The  points  are  ordered  as  indicated  in  Table  H# 


Table  II 


Sequence  of  points 


0  0 

1 
2 


pz2  ♦  1 
0 


pa2  ♦  1 


• 


• 


pzl+1  pz2+l 

As  explained  in  the  Procedure  the  probability  of  scattering 
from  group  to  group  and  fission  cross  section  are  associated  with  the 
grid  points*  These  are  also  placed  on  the  tape  of  data  immediately 
following  the  flux  for  f  *  3,  i  ■  2,  and  /*  2  in  the  order 
indicated  on  the  following  page  in  Table  HI. 
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Table  HI 
Sequence  of  additional  data  for  a  point 

T21 

r31 

22 

r 

r32 
r33 

The  total  scattering  cross  section  is  associated  with  the  region 
between  points  and  are  ordered  as  in  Table  ]V# 

Table  17 
Sequence  of  Total  cross  sections 

1  H 

0  0 

1 

• 

2 

pzl  0 

pz2 
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2 

fc  TAPE  I4IO-3I6-IO32  S     AND  GAUSS  CAMPBELL  AND  PANCIERA 

(22|,6) 

pzl-5     pz2«5  pz3=20    pzi^-20    pfl«5000  pf2*10000  pf3»61*li 

pflt»36  pf 5-701*  pflO-2650  pfl3-lli70  pfl£-300 

/Constants 

al.,+1.0        -0.333333  +O.6666667  a2,+0.0      +1.333333  +2.0 

a3,+2.1j6        *U,-1.0  a£,+1.96770#i  +1.1306938 

♦1.306938      +1.967705U  a6,+0.03    +0.0?  +10.0  a7,*0.3l83 

+0.3631*  +0.3183  a8,+0.652ll£2  +0.3li78&8  *0.3ltf85U8 

♦0.65211*52    a9,+22.0  alO,+U.O    all, +1.691*1*625  +0.381*331*1 

/Cross  Sections 

bl5*1.0         DITTO  12r/b2,+1.0  DITTO  12r/m6,+0  m7,+0 

m8i)pf3  m9,pfl*  mlO,+92       mll,+6        ml2,pfl      ml3,pf2      mll*,+0 

ml5*+0 
/RADII 

r3j>+0.0  rl*,+0.0  r5,+0.0 

/NORMALIZING  CONSTANTS 

dl>+0.0  d29«K).0  d3,+0.0      rl,+0.0      r2,+0.0       slO,+0.0 

Aariables 

el.,+l*.8   c25-5.6  c3,+0.6885        +0.2295   +0.1*590   +3«0355 

♦1.0118    +2.0236 

/j°lj,m-l  numbers  m 

nl^+1.0   DITTO  72r/sl,+1.0       +1.0     +1.0     fl,+1.0 

♦1.0       +1.0  +1.0     +1.0     +1.0     +1.0 

/j^l^m  numbers 

n29+1.0  DITTO  72r/s2,+1.0  .  +1.0  +1.0  t2,+1.0 

+1.0  +1.0  +1.0  +1.0  +1.0  +1.0 

/j^l^m+l^,  numbers 

n35+1.0    DITTO  72r/s3,+l«0       +1.0     +1.0     t3,+1.0 

♦1.0       +1.0  +1.0     +1.0     +1.0     +1.0 
/zs   m-l^,  numbers 

nl^+1.0    DITTO  72r/sl*,+1.0        +1.0     +1.0     tl*,+1.0 

+1.0       +1.0  +1.0     +1.0     +1.0     +1.0 

/^m  numbers 

n5<>+1.0  DITTO  72r/s5,+1.0  +1.0  +1.0  t5,+1.0 

♦1.0  +1.0  +1.0  +1.0  +1.0  +1.0 

/jgm+lj,  numbers 

n65+1.0  DITTO  72r/s6,+1.0  +1.0  +1.0  t6,+1.0 

+3U0  +1.0  +1.0  +1.0  +1.0  +1.0 
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/j+l5m-l,  numbers 


n75+loO 

DlTTTO 

72r/s7, 

+1.0 

+1.0 

+1.0 

t7,+1.0 

+1.0 

+1.0 

♦1.0 

+3,.0 

+1.0 

+1.0 

/j+l.ra,,  numbers 

n85+l«0 

DITTO 

72r/s8, 

+1.0 

+1.0 

+1.0 

t8,+1.0 

+1.0 

+1.0 

+1.0 

+1.0 

+1.0 

+1.0 

/j+l^m+l,  numbers 

n9*+1.0 

DITTO 

72r/s9, 

+1.0 

+1.0 

+1.0 

t9,+1.0 

+1<>0 

+lo0 

+1.0 

+1.0 

+1.0 

+1.0 

/CONVERGED  GAMMA 

ml<,  itaml8 

itsrS 

idva9 

itscl 

iscO 

icr3 

iscl 

icr6 

isc2 

icr2 

m2jicac2 

idvcl 

iac2 

idvall+c 

iscO 

imral+c 

icsl 

itsc3+c 

iscO 

ictm2 

isc2 

ictra2 

iscO 

icrpz3 

m3,icaal 

itsslO 

icacl 

itsrl 

Vise! 

icrpzl+1 

isc2 

icrpz2+l 

mii^ispil 

isc2 

ictmli 

icarl 

itsr2 

iadcl 

iadcl 

itsrl 

iscl 

ictmli 

isep 

ietn£ 

iSTOP 

m^icaslC 

•  icpm3 

ml8^ispo 

/3,m  SENSE  AND  CONTROL 

il,,itai7 

iscl 

itim6 

isc2 

itim7 

OUT 

cam? 

mhralO 

slhl£ 

adml2 

tsmlS> 

cam8 

mhm6 

slhl5 

tsmlii 

adml5> 

tsi6 

caml£ 

sumlO 

tsml£ 

admlli 

tsm29 

camlli 

sum8 

adnk5> 

tsm26 

camlli 

adm8 

adml£ 

tsm30 

cam7 

sul 

mhmll 

slhlS 

aflml3 

tsml5 

cam9 

mhmo 

slhl5 

tsmlii 

admlf> 

tsm32 

camlli 

sum9 

adml5 

tsm31 

cam6 

dmO 

cpi2 

cam? 

dmO 

cpi3 

IN 

i»pm23 

ispm27 

ispi5> 

12 9 cam? 

dmO 

epili 

IN 

ispra23 

ispm21 

ispif> 

135IN 

ispm23 

ispm2U 

ispi£ 

±k,TS 

ispm23 

ispml6 

i^iDOB 

n$ 

i6,+0 

♦78 

i?jispo 

/DATA  READ 

IN  PROGRAMME 

m^itamSa 

iDIB 

nl 

ra26,+0 

+276  ( 

iDIB 

nii 

m29s+o 

+2?6 

iDIB 

n7 

m30,+O 

♦276 

iDIB 

bl 

m315+0 

+6 

iDIB 

b2 

m32,+0 

♦6 

m332ispo 

/ORIGIN  PROGRAM 

Ifil^itaral? 

ispplO 

isppl3 

icaal 

idvt7 

itsdl 

isp«3 

iDIB 

pfl3 

+20^0 

♦17S 

ispfl 

ml7j>ispO 
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isppl3 
iDIB 


isppl3 
iDIB 


/CENTERLINE  PROGRAM 

m21aitam22  isppli 

+175       ispfl 

m225lsp0 

/MIDPLANE  PROGRAM 

m2ii3itam25  ispp7 

+175       ispfl 

m255ispQ 

/GENERAL  POINT  PROGRAM 

m275itaiB28  isppl    isppl3 

+175       ispfl    iDIB 

iDIB       pfl3     +2350 

+21*75      +175     ispfl8 

/Source  S°5S15S2  at  all  Point 

pl35itaplh  icat6    imrt5 

iadtl      itst7    icat6 

iadtl      itst7*2  j  -icat6 

iadtL=.     itstl    icat6+l* 

/Source  Normalization+Tolerance  Program 

•3sitae7    icat7+2   irardl    itss5+2 

ieat7  inrrdl  itst7 

iepe5  e69iesa6    itsslO 

e5pieat7        itss5  e7<,ispO 

/SCALAR  Flux-General  Point 


ispe3 
pfl3 


ispe3 
pfl3 


ispe3 
pfl3 
+125 
ra28,ispO 

itstl 

imrt5+2 
imrt5"fli 
imrt5+10 


isus5 
ispe5 


iDIB 
+2ltf5 


iDIB 
+2225 


iDIB 

+2225 

ispf27 


icaa3 
itstl 
itstl 
iadtl 

icat7+l* 

icpeii 

el^idvt7 


icaa7 
iadt2 
itstl 


pl5itap3  isc3     icr3 

imrn5+6+c  iactl  itstl 

imr&8  itstl    icaa7 

imrn5+2l*+c  iadt2  itst2 

imra8+2  iadtl    itstl 

ieaa7*2  imrn5+l|2+c 

iadt2  imra8+li      iadtl 

itst2  icaa7+2      iram5+60+c 

iram5+66+e  iadt2  imra6+6      iadtl 

/Scalar  Flux°Centerline  Point 

pi*5itap6  isc3  icr3  p5picaa8 

imm5+l8+c  iadtl  itstl  icaa8+l* 

itstl  icaa8+6      imrn5+5li+c 

ieaa8  imrn5+6+c  itst2    icaa8+2 


pfl3 
♦175 


pfl3 
+125 


pfl3 
+125 

iDIB 


imrt5*12 
icat6+2 
icat6+2 
itst7+U 

imrdl 
idvt7 
iada6 


itstl 


+2050 
ispfl8 


+2050 
ispflO 


+2050 

ispflO 

pfl3 


imrt6+lt 
imrt5+6 
imrt5+8 
plli9isp0 

its  s  5+1* 

isua6 

icpe6 


p2sicaa7  imrn5+c 

icaa7+ii      imrn5+12+c 

imrn5+l8+c  itst2 

icaa7+li  irarn5+30+c 
imrn5'f36+c 
itst2  icaa7+l* 

icaa7  imrn5+51i'fc 


iact2 
itst6+c 


itst2 
ictp2 


imrn5+c   itstl 
imrn5+36+c 
iadtl    lmra7 
imrn5+2l*+c 


icaa7+2 

iadtl 

icaa7+2 

iadt2 

itst2 

imrn5+lt8+c 

icaa7+l* 
P3,i»p0 

icaa8+2 
iadtl 
itstl 
iadt2 
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itst2              icaa8*ii 

imm£+li2+c 

iadt2 

itst2 

icaa8+6 

imrn5+60+e     iadt2 

imra7*2 

iadtl 

imra2+l| 

itst6+c 

ictp5 

p65isp0 

/Scalar  Flux-Mi dplane 

i  Point 

p7^itap9        isc3 

icr3 

P8,icaa7 

imrnS^c 

itstl 

icaa7+2 

imrn5*6*c       iadtl 

itstl 

icaa7*U 

irarn£+12+c 

iadtl 

imra8              itstl 

icaa7 

imm5+l8-»c 

itst2 

icaa7*2 

£ar&£+2b+c    iadt2 

itst2 

icaa7+li 

imrn!?+30+c 

iadt2 

imra8+2          iadtl 

imra2+ii 

itst6+c 

ictp8 

p9,ispO 

Scalar  Flux  -  Origin 

plO^itap^    isc3 

icr3 

pll,icaa7  inarnS>+c 

itstl 

icaa7*2 

imrn^+o^c       iadtl 

irara8 

itstl 

icaa7 

iinrn^lS+c 

itst2              i©aa7+2 

imrn£+2lr»-c 

iadt2 

imra8«-2 

iadtl 

irara2+lt          imra2*]i 

itst6+c 

ictpll 

P12,ispO 

/Radius  extrapolation  routine 

qij>itaq9        icad2 

isual 

itstl 

icpq2 

q3,itst3 

icad3 

isual              itst2 

icpqli 

q5,  itstl* 

isut3 

icpq6 

icad2 

i»xd3              itsd2 

icar3 

iexr/i 

itsr3 

icat3 

iextU 

itst3              icatl 

iext2 

itstl 

q6,icat3 

isua6*2 

icpq7 

q8,icad2         isud3 

itst6 

icarli 

isur3 

imrtl 

idvt6 

iadr3              itsr3 

q9,isp0 

q2,imral| 

ispq3 

q^irarali 

ispqf> 

/Control 

xl5i©arli        iscli 

icrpzli 

ispml 

icadl 

itsd3 

x2,icar3 

ispml              icadl 

itsd2 

iscli 

ictx3 

iSTOP 

x3,ispql 

ispx2               q7,iDIB 

Pf5 

pflO 

pfl5 

qlO,isppf£ 

iSTOP 

/Forking  flux 

lltf0/yl5+0    DITTO 

17^r/y2, 

+0 

/First  Quadrant  Flux 

/First  Flux 

DA20£0/  litfO/  fl,itaf3 

isc3 

icr3 

f2,icsb2*6-»'C 

imra5>              iaircl 

itstl 

iadal 

imrn9+c 

itst2 

icas5'*,c 

iads9+e          imra^ 

irarcl 

iadt2 

irarc3 

itst2 

icstl 

irar^3              iadal 

itstl 

icaal 

isuc3 

imrn6*c 

iadt2 

id^fel              itsn£+c 

/Second  Flux 

fij.giesa2+2     idvrl 

itstl 

icsb2+6+c  imra£ 

itst2 

iadtl 
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imrcl 

itst3 

iadal+2 

imrn9+6+c  itstl; 

icitl 

iadt2 

imrcl 

itst6 

iadal+2 

imrn9+c   iadtl; 

itstl; 

icat6 

isual+2 

f59imrn5>+c 

iadtl;    itstl; 

icas5+c 

iad»9+c 

imra^ 

imrcl 

imra2+l; 

iadtl;    itstl; 

icsal+1; 

iadal+2 

id"*a2+l; 

itstl 

ican6+c 

iadn£+c   isun9+c 

isun8+c 

imrtl 

iadtl; 

imrc3+2 

itstli 

f6,icst3  imrc3+2 

iadal+2 

itst2 

icaal 

isuc3+2 

iuiral+2 

itst3    ican6+6+c  iadn6+c 

isun£+c 

imrt3 

iadtl; 

idvt2 

itsn^+6+c 

/Third  Flux 

f7«,iesb2+6+c 

irara5+2 

imrcl    itstl 

iadal 

imrn9+l8+c 

itst2 

icas£+c 

iads9+c 

imra5>+2   imrcl 

iadt2 

irarc3+6 

itst2 

icstl 

iadal 

imrc3+6   itstl 

icac3+6 

isual 

itst3 

ican8+l8+c 
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